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■ Abstract. Spectroscopy in the infrared provides a means to assess im- 
| portant properties of the plasma in gaseous nebulae. We present some 

_ ■ of our own work that illustrates the need for interactions between the 

', themes of this conference - astronomical data, atomic data, and plasma 

jy^ • simulations. We undertook Infrared Space Observatory (ISO) observa- 

, tions with the intent of better understanding the effects of density vari- 

■ ations in nebulae, particularly planetary nebulae (PNs), by determining 
average electron densities from the flux ratios of several fine-structure, 
IR emission lines. Instead, we are able to ascertain only minor density 

, information because of several instances of the observed line flux ratios 

00 

; being out of range of the theoretical predictions using current atomic 

ff) . data. In these cases, the ISO data cannot presently be used to derive 

electron density, but rather provide direction for needed improvements in 
the atomic collision strengths. 
^ I We have detected an unidentified (uid) strong emission line in an 

ISO/SWS spectrum of the Orion Nebula. The line has a rest wavelength 
2.89350±0.00003 /im. A long-slit UKIRT observation confirms the pres- 
ence of this line and shows that the emission is spatially extended and 
appears to be coincident with the brightest part of the ionized region. 

■ We do not detect the uid line in our SWS02 spectra of any of the several 
| bright PNs which we observed for a comparable time. The need for ba- 
sic atomic data, in this case wavelengths to aid species identification, is 

• i-H ! paramount for future progress. 

^ ' We look toward the future with a brief synopsis of upcoming or 

^ . planned IR missions that promise a significant spectroscopic bearing. The 
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missions (and their instruments) discussed are SIRTF, SOFIA, and HER- 
SCHEL (upcoming) and NGST, GSMT, and SAFIR (planned). 



1. Introduction 

Most observational tests of the chemical evolution of the universe rest on emis- 
sion line objects; these define the endpoints of stellar evolution and probe the 
current state of the interstellar medium. Gaseous nebulae are laboratories for 
understanding physical processes in all emission-line sources, and probes for 
stellar, galactic, and primordial nucleosynthesis. 

There is a fundamental issue that continues to be problematic - the discrep- 
ancy between heavy element abundances inferred from emission lines that are 
collisionally excited compared with those due to recombination/cascading, the 
so-called "recombination lines". Studies of planetary nebulae (PNs) contrast- 
ing recombination and collisional abundances (Liu et al. 1995, Kwitter & Henry 
1998) often find differences exceeding a factor of two. In an extensive study 
of NGC 6153, Liu et al. (2000) found that C++/H+, N++/H+, 0++/H+, and 
Ne ++ /H + ratios derived from optical recombination lines are all a factor of ~10 
higher than the corresponding values deduced from collisionally-excited lines. 
Abundances determined from these two methods disagree by a factor larger 
than the spread of abundances used to determine such fundamental quantities 
as Galactic abundance gradients (e.g., Shaver et al. 1983; Simpson et al. 1995; 
Henry & Worthey 1999). 

Most of the efforts to explain the abundance puzzle between collisional and 
recombination values have attempted to do so by examining electron tempera- 
ture (T e ) variations in the plasma. This is often done, using the formalism of 
Peimbert (1967), in terms of the mean-square variation (t 2 ) of T e . The inferred 
metallicity obtained by using the usual (optical/UV) forbidden lines is very sen- 
sitive to T e (exponential) and t 2 . On the other hand, recombination lines are 
rather insensitive to T e and t 2 . Agreement close to the higher recombination 
value can be forced in the derived abundance by attributing the difference to 
(solving for) t 2 . For instance, consider the case of the PN NGC 7009, which has 
stood near the center of the abundances controversy. Liu et al. (1995) found 
that the recombination C, N, and O abundances are a factor of ~5 larger than 
the corresponding collisional abundances. By invoking t 2 ~ 0.1 for oxygen in 
NGC 7009, agreement can be forced in the derived abundance close to the higher 
recombination value - a value more than 2.5 times larger than the solar O/H of 
7.41xl0~ 4 (Grevesse Sz Sauval 1998). Such a large t 2 is not at all predicted by 
current theory/models (e.g., Kingdon k, Ferland 1998). 

For NGC 7009, Rubin et al. (2000) produced a T e map on a pixel scale of 
0.1" from the flux ratio of [O in] 4363/5007 using a set of HST/WFPC2 images 
that included narrow-band filters F437N and F502N to assess these line fluxes. 
Their analysis of t 2 in the plane of the sky indicates very small values, <^ 0.01, 
throughout the nebula. It seems clear that something more than T e variations 
alone is necessary to explain the abundance dilemma. 

We have been addressing other ingredients that may affect the determina- 
tion of elemental abundances from observations. Perhaps the most important 
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Figure 1. Three valuable indicators of N e - the intensity ratio of 
[N n] (205/122), [O in] (88/52), and [S III] (33/19) vs. N e (cm- 3 ). The 
[N II] ratio is most sensitive at the lowest N e , and [S ill] at the highest. 

of these is density variations. While T e variations are expected theoretically to 
be fairly small in PNs and H II regions, variations in electron density (N e ) are 
almost certainly very large. High N e values (~10 6 cm -3 ) are inferred in the 
cometary knots of the Helix nebula (e.g., Walsh & Meaburn 1993; Burkert & 
O'Dell 1998) and around proplyds in Orion (e.g., Henney & O'Dell 1999 and 
references therein). Gas at N e = 10 6 cm -3 will produce 10 6 times as much 
emission as the same volume of "normal" nebular gas at N e = 10 3 cm -3 for re- 
combination lines and collisionally-excited lines that are not suffering collisional 
deexcitation - when N e < N cr u (the critical density). 

Several well known line ratios of collisionally-excited lines serve as diagnos- 
tics of N e , probing different ionization conditions and different density regimes. 
Because various collisionally-excited lines undergo collisional deexcitation at dif- 
ferent N cr it values, the derivation of an average N e for a given line ratio depends 
on the specific N e (and T e ) dependence of the respective volume emissivities. 
In general, different average N e values are obtained for various line pairs due 
to different contributions to the observed intensities from the volume observed. 
Very substantial biases in the inference of abundances from collisionally-excited 
lines may then result (Rubin 1989; Viegas & Clegg 1994; Liu et al. 2000). 

2. The Promise of the Infrared for Density Diagnostics 

Measurement of the two AJ = 1 IR fine-structure transitions for species with a 
3 _p ground state provides a diagnostic for N e . The ratio of their fluxes is a sen- 
sitive indicator of N e (over some range in N e ) that generally is not significantly 
affected by T e or extinction. For example, the flux ratios [O in] 52/88 fim, [S in] 
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Figure 2. Similar to Fig. 1 for two other indicators of N e - the inten- 
sity ratio of [Ne m] (36/15) and [Ar m] (22/9). These diagnostics are 
now available from ISO. These ratios are sensitive indicators of higher 
density material than those ratios in Fig. 1. 

19/33 /mi, and [N II] 122/205 /mi provide N e [0 III], N e [S III], and iV e [N II], 
respectively, as shown in Figure 1 (Rubin et al. 1994). 

There is only a weak dependence on T e as illustrated for each species by the 
three curves for T e = 5000, 8000, and 10000 K. For conditions prevalent in many 
PNs and H n regions, these lines readily show the effect of collisional deexcitation 
because some of the plasma may exceed their respective N cr n- values. We are 
able to use only two of the above diagnostics because ISO did not observe the 
205 /im line. 

With ISO, there are now at least three important additional IR line ratios 
available - [Ne III] 36.0/15.5, [Ne v] 14.3/24.3, and [Ar III] 21.8/8.99 /mi - for 
studies of N e structure and variations. This is because ISO was able to observe 
the [Ne in] 15.5 /mi, [Ne v] 14.3 /mi, and [Ar in] 21.8 /mi lines, which suffer very 
seriously from the atmosphere even for airborne astronomy, e.g., for the Kuiper 
Airborne Observatory (KAO) or SOFIA. (Observations of these lines with the 
IRAS LRS were sparse and of dubious quality.) These three ratios are sensitive 
indicators of higher density material than are the other three ratios mentioned 
earlier. Two of these ratios - [Ne in] 36/15 /mi and [Ar in] 22/9 /mi - are 
displayed in Figure 2. 

The ratios are similar in that they are most discriminant in the range 3.8 <J 
log N e <^ 5.4. We will use two more iV e -diagnostic line ratios [Mg v] 13.5/5.61 /mi 
and [Ar v] 13.1/7.9 /im. These will be used to probe the higher ionization gas. 
Although usually weaker than the lines mentioned earlier, all are expected to 
be seen in a large number of high ionization PNs. The hope for the ISO IR 
data is that these numerous "new" N e diagnostics, in combination with the 
N e diagnostics available from the optical/UV, will permit a viable tomographic 
analysis of the N e structure in PNs. 
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3. ISO Observations and Results for the Density Analysis 

Under our GO programs, we observed the PNs NGC 2022, NGC 6210, NGC 6818, 
and IC 2165 with the ISO Short Wavelength Spectrometer (SWS). Our obser- 
vations with SWS were made in SWS02 mode, which provides higher spectral 
resolution than SWS01 mode. One of our goals was to derive N e from several di- 
agnostic line pairs in order to address density variations within these objects. As 
described above, we aimed to obtain N e values from the following six flux ratios: 
[S in] (18.7/33.5), [Ar III] (8.99/21.8), [Ne III] (15.5/36.0), [Ar v] (7.90/13.1), 
[Ne v] (14.3/24.3), and [Mg v] (5.61/13.5/mi). Because the ISO/SWS aperture 
size depends on wavelength, the above line fluxes are not always directly compa- 
rable. For the lines discussed here, the sizes were: 14" x 20" for the 5.61, 7.90, 
& 8.99 fim lines; 14" x 27" for the 13.1 - 24.3 fim lines; and 20" x 33" for the 
33.5 & 36.0 /mi lines. Table 1 provides a summary of which lines were observed 
in the various sources. 



Table 1. 



ID Vac. A (imi) 

6 5.60985 

4 7.90160 

2 8.99138 



4 13.1022 

6 13.5213 

5 14.3217 

3 15.5551 



2 21.8293 
5 24.3175 



Species (terms) 

Aperture 14" 
Mgv ( 3 P t - 3 P,) 

Arv ( 3 Pj - 3 Pi) 

Arm ('Pr - 3 P 2 ) 



Source SWS band 
X 20" 

102165 2B 
NGC2022 

IC2165 2C 
NGC6818 

IC2165 2C 
NGC2022 
XGC6210 
NGC6818 



Aperture 14" 

Ar v ( 3 P t - 3 P ) 

Mgv ( 3 P„- 3 Pi) 
Nev ( 3 P 2 - ^Pr) 

Nem f'Pr - 3 Pj) 

S III ( 3 Pj - 3 Pi) 

Arm ( 3 P„ - 3 P t ) 
Nev ( 3 Pr - :1 P ) 



X 27" 

IC2165 

NGC2022 

NGC6818 

IC2165 
NGC6818 

IC2165 

NGC2022 

NGC6818 

IG2165 
NGC2022 
XGC6210 
NGC6818 

IC2165 
XGC2022 
XGC6210 
NGC6818 

NGC6210 
NGC6818 

IG2165 
XGC2022 
XGCG210 
NGC6818 



3 36.0135 



Aperture 20" 

S III ( 3 P t - 3 P ) 



Nem ( J P - J Pi) 



X 33" 

IC2165 
NGC2022 
XGCC210 
NGC6818 

IG2165 

NGC6210 

XGCC818 



The position angles (PA) of the long-axis of the various apertures were 
183.17, 13.27, 349.07, and 193.71° for NGC 2022, NGC 6210, NGC 6818, and 
IC 2165, respectively, and rough estimates of the respective source "diameters" 
are 19", 16.2", 20", and 9". However, there is certainly emission well beyond 
these "diameters" in the case of NGC 6 210 see an overlay 
( http:/ /www-space.arc.nasa.gov/^rubin/turtle.html ) of our ISO apertures on 
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Figure 3. ISO SWS02 spectrum of [Ne v] 14.32 /mi line in the PN 
NGC 6818. This important ^-diagnostic line cannot be observed even 
from an airborne platform. The data (asterisks) have been fit with a 
Gaussian and a linear baseline. The residuals, data point minus curve 
fit, are indicated as the dashed line. 

HST/WFPC2 images taken by us under GO-6 792 - and NGC 6818 
dhttp: / / www-space.arc.nasa.gov/^rubin / zorro.html ) other WFPC2 images un- 
der GO-6792. The best object for our analysis in this sense is IC 2165, which is 
smaller than the smallest aperture. For the other PNs, the use of line flux ratios 
may still be valid as long as the specific ionic emitting zone is circumscribed by 
the smaller of the two apertures. Of the six line sets, only the [Ne v] pair of 
lines (14.3, 24.3) was observed with the same aperture size. 

Figures 3 and 4 are ISO SWS02 spectra of the [Ne v] 14.32 and 24.32 /mi 
lines in NGC 6818. The aperture used for each was 14" x 27". The data are fit 
very well by a Gaussian profile. Data processing was performed using the ISO 
Spectral Analysis Package (ISAP). NGC 6818 is not fully enclosed by a 14" x 27" 
aperture. However, the same portion of the nebula is observed in both of these 
lines. Because the high ionization Ne +4 zone is much more concentrated toward 
the central star than the entire ionized (H n) region is, it is highly probable that 
the total Ne +4 zone is enclosed within the observed aperture. 

Figure 5 shows the theoretical flux ratio F(14.3)/F(24.3) vs. N e (cm -3 ), 
using effective collision strengths from Lennon & Burke (1994). Throughout 
this paper, we shorten effective collision strength (T) to collision strength (CS). 
Our ISO data with regard to this ionic species are most interesting. In the 
subsections to follow, we present the SWS data and discuss the six ionic species 
in order of lowest to highest degree of ionization as indicated by the range in 
ionization potential (LP.) over which the species exist. 



3.1. S m 

Observations were made of both the [S in] 18.7 and 33.5 /mi lines in all four of 
our program PNs. The line profiles were fit with the line-fitting routine in ISAP, 
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Figure 4. Similar to Fig. 3 for the ISO SWS02 spectrum of [Ne v] 
24.32 jum line in the PN NGC 6818. 
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Figure 5. An outstanding indicator of N e is the flux ratio of [Ne v] 
(14.3/24.3). This shows the theoretical flux ratio F(14.3)/F(24.3) vs. 
N e (cm -3 ) using Lennon Sz Burke (1994) collision strengths. The ratio 
is insensitive to T e (10,000 K solid, 15,000 K dashed curve). This 
diagnostic is now available from ISO data. 
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which provides the line flux and uncertainty. We fit the data points with a 
Gaussian and a linear baseline such as shown in Figs. 3 and 4 for the [Ne v] 
data. In general, the data were fit well by this method. Table 2 presents the 
measured fluxes and uncertainties. These uncertainties do not include systematic 
errors. The last column has the N e derived from the flux ratio assuming a T e 
of 10,000 K. It is likely that only the N e for IC 2165 is valid. Because of the 
extended [S m] emission and the mismatch of aperture sizes, the F(18.7)/F(33.5) 
ratio is likely a lower limit and the N e entry a lower limit for the other objects. 

In the lower part of Table 2, we provide the theoretical F(18.7)/F(33.5) 
ratio in the low-iV e limit; there, the smallest value of the ratio obtains. The 
transitions are shown in a term diagram (not to scale) with 7V cr .j t -values for 
levels 2 and 3 given for T e = 10,000 K. 

[S III] - I.P. range 23.33 - 34.83 eV 



F(18.7) (W cm- 2 ) Ff33.5) F(18.7)/F(33.5) N, (cm" 3 )* 



Table 2. 



Source 
IC2165 

NGC2022 

NGC6210 

NGC6818 

* for T c = 10,000 K. 

*********************************************************** 

Theoretical Ratio in Low-A r c Limit: lowest value at low-A r e limit 

Ngrit (cm- 3 ) 



4.879E-19 
±7.656E-21 



1.561E-19 
±8.529E-21 



2.196E-18 
±3.221E-20 



1.475E-18 
±2.014E-20 



2.048 E-19 
±1.747E-20 



1.659E-19 
±8.633E-21 



1.111E-18 
±3.473E-20 



1.326E-18 
±3.181E-20 



2.382 
0.941 
1.977 
1.112 



3220 
539 

2400 
823 



% IK) Ratio 

10,000 0.6030 

15,000 0.6757 

20,000 0.7328 





18.7 fim 


33.5 /j,m 





3.2. Ar in 

Observations were made of the [Ar m] 8.99 and 21.8 /im lines in NGC 6210 and 
NGC 6818. Table 3 presents the measured fluxes for these data. The comments 
regarding line fitting, etc., and the tabular entries that accompanied Table 2 
apply here and for the ionic species to follow. In the lower part of the Table, 
we have the theoretical F(21.8)/F(8.99) ratio in the low-iV e limit, where the 
highest value of the ratio obtains. For NGC 6210, the observed ratio is close 
to the low-Af e asymptotic limit, where the ratio is insensitive to N e and thus 
not useful for deriving N e . For this and other such instances, "near low limit" 
is entered in the A^-column of the Table. For NGC 6818, the observed ratio 
exceeds the low-AT e limit for any reasonable T e . Because the [Ar III] emission is 
likely more extended than the smaller 14" x 20" aperture used for the 8.99 /im 
line, the F(21.8)/F(8.99) ratio is likely an upper limit and thus any inference 
that the observed ratio is out of the theoretical bounds is uncertain. We enter 
"XXX low ?" in the Table. 
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Table 3. 



Source 
NGC6210 



[Ar III] - I.P. range 27.63 - 40.74 eV 

F(8.99) (W cm- 2 ) F(21.8) F(21.8)/F(8.99) N„ (erg- 1 ) 

0.0670 near low limit 



9.676E-19 
±1.754E-20 



5.465E-19 
±1.349 E-20 



6.478E-20 
±6.777E-21 



4.806E-20 
±6.682E-21 



*********************************************************** 

Theoretical Ratio in Low-A r e Limit: highest value at low-A^ limit 

N m -, (tmi 



T„ (K) Ratio 

10,000 0.0680 

15,000 0.0704 

20,000 0.0720 





21.8 fim 


8.99 fim 


3 P, 



3.3. Ne in 

Observations were made of the [Ne m] 15.5 and 36.0 /im lines in IC 2165, 
NGC 6210, and NGC 6818. Table 4 presents the measured fluxes for these 
data as well as the theoretical F(36.0)/F(15.5) ratio in the low--/V e limit, where 
the highest value of the ratio obtains. For all the sources, the observed ratio 
exceeds the low-iV e limit. However because the [Ne III] emission for NGC 6210 
and NGC 6818 may be more extended than the smaller 14" x 27" aperture used 
for the 15.5 /im line, the F(36.0)/F(15.5) ratio is likely an upper limit and again 
has an uncertain (XXX low ?) conclusion. On the other hand, for IC 2165, 
we are dealing with integrated fluxes and conclude robustly that the observed 
ratio is out of the theoretical bounds and enter "XXX low" in the Table. The 
theoretical ratio in the low-iV e limit is set by the effective collision strengths. 
We suggest the data for IC 2165 are pointing to a need to reexamine the CS, 
particularly those between the 3 P fine-structure levels. 



Table 4. 



Source 
IC2165 

NGC6210 

NGC6818 



[Ne III] - I.P. range 40.96 - 63.45 eV 
F(15.5) (W cm' 2 ) F(36.0) F(36.0)/F(15.5) N n |im-'| 



3.130E-18 
±3.031E-20 



1.614E-17 
±2.700E-19 



4.852E-18 
±3.896E-20 



3.097E-19 
±6.295E-21 



1.666E-18 
±6.222E-20 



7.057E-19 
±1.629E-20 



0.0989 



XXX low 
XXX low ? 
XXX low ? 





36.0 /im 


15.5 /j,m 


3 P, 



y„n (cm 3 ) 
2.92E4 



*********************************************************** 
Theoretical Ratio in Low-A r e Limit: highest value at low-A r e limit 
T r (K) Ratio 
10,000 0.0879 
15,000 0.0873 
20,000 0.0860 
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3.4. Ar v 

Observations were made of the [Ar v] 7.90 and 13.1 lines in IC 2165 and 
NGC 6818. Table 5 presents the measured fluxes for these data as well as the 
theoretical F(7.90)/F(13.1) ratio in the low-iV e limit, where the lowest value of 
the ratio obtains. For both PNs, the observed ratio is near the low-A^ limit. 
Here, even in the case of NGC 6818, it is likely that the Ar +4 zone is much 
more centrally concentrated than the size of the H + zone; but if the [Ar v] 
7.90 emission were more extended than the smaller 14" x 20" used, then the 
F(7.90)/F(13.1) ratio would be a lower limit. 



Table 5. 



Source 
IC2165 



[Ar v] - I.P. range 59.81 - 75.04 eV 

F(7.90) (W cm- 2 ) F(13.1) F(7.90yF(13.H JV„ (cm-'l 

0.611 near low limit 



5.954E-20 
±3.904E-21 



1.057E-19 
±1.891E-21 



9.748E-20 
±1.844E-20 



1.833E-19 
±4.826E-21 



near low limit 



*********************************************************** 
Theoretical Ratio in Low-A r e Limit: lowest value at low-A r e limit 

N„«. (cm- 3 ) 
1.60E5 



T r m 


Ratio 








10,000 


0.5724 








15,000 


0.6134 








20,000 


0.6411 




7.90 iivn 






13.1 (jm 





p? 



Table 6. 



Source 
IC2165 



[Ne v] - I.P. range 97.11 - 126.21 eV 
F(14.3) (W cm- ') Ff24.3) F(14.3)/F(24.3) jV, (cm' 3 ) 



2.284E-18 
±3.838E-20 



1.957E-18 
±1.844E-20 



* for T. = 10,000; 1950 cm 3 at 15,000 K. 



NGC2022 1.319E-18 
±1.982E-20 



NGC6818 3.701E-18 
±3.384E-20 



1.596E-18 
±2.108E-20 



5.233E-18 
±5.306E-20 



0.827 
0.707 



XXX low 
XXX low 



Theoretical Ratio in Low-A r c Limit: lowest value at low-A r e limit 

Neril (CHI-') 



T,. IK) Ratio 

10,000 0.9067 

15,000 0.8935 

20,000 0.8814 

25,000 0.8712 





14.3 ,11m 


24.3 fim 





3.74E4 



3.5. Ne v 

Observations were made of the [Ne v] 14.3 and 24.3 fim lines in IC 2165, 
NGC 2022, and NGC 6818 (see Figs. 3 k. 4). Table 6 has the measured fluxes for 
these data as well as the theoretical F(14.3)/F(24.3) ratio in the \ow-N e limit, 
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where the lowest value of the ratio obtains. For NGC 6818, the observed ratio is 
well below the low-iV e theoretical limit for any reasonable T e ; for NGC 2022, the 
observed ratio is also out of bounds. For IC 2165, we derive N e = 1310 (if T e = 
10,000 K) and 1950 cm' 3 (if T e = 15,000 K) from Fig. 5 (see Table 6). Com- 
pared to iV e -values derived for IC 2165 by Stanghellini & Kaler (1989), albeit 
for other ions, the values here are very low. 

We are aware of three other measurements of N e from the [Ne v] ratio 
using ISO data. For the Seyfert galaxy NGC 4151, Alexander et al. (1999) 
found an observed ratio that is close to the low-iV e limit For the PNs NGC 7027 
and NGC 6302, van Hoof et al. (2000), found N e of 26,900 and 12,300 cm~ 3 
using T e ~20,000 K. Compared to A^ e -values derived from optical data for these 
PNs for other ions (Stanghellini &: Kaler 1989), the IR values are lower. The 
fact that these previously published [Ne v] ratios and our new IC 2165 result 
yield generally lower iV e -values than other methods tends to lend support to our 
conclusion that the low-iV e theoretical limit may need re-examination per our 
XXX low results. Just because the NGC 7027, NGC 6302, and IC 2165 data 
lead to "legal, in-theoretical-bounds values" does not mean the inferred N e s are 
correct and does not rule out the possibility that the collision strengths between 
the 3 P levels may need revision. 

3.6. Mg v 

Observations were made of the [Mg v] 5.61 and 13.5 u,m lines in IC 2165. Table 7 
presents the measured fluxes as well as the theoretical F(13.5)/F(5.61) ratio 
in the low-N e limit, where the highest value of the ratio obtains. Again, the 
observed ratio is out of the theoretical bounds. 



Table 7. 



Source 
IC2165 



[Mgv] - I.P. range 109.31 - 141.27 eV 

F(5.6H (W cm' 2 ) F(13.5) F(13.5)/F(5.61) N„ (cm' 3 ) 
0.1032 XXX low 



2.692E-19 
±9.792E-21 



2.779E-20 
±4.807E-21 



*********************************************************** 
Theoretical Ratio in Low-A r e Limit: highest value at low-A r a limit 



T„ 0O Ratio 

10,000 0.0797 

15,000 0.0814 

20,000 0.0819 



3 P 





13.5 fim 

3 p, 


5.61 /im 


3 p, 



4. Observations of a Strong Unidentified IR Emission Line 

We have detected an unidentified (uid) strong emission line in an ISO SWS02 
spectrum of the Orion Nebula that was obtained for our Guest Investigator 
program on 1997 October 12 (UT), see Figure 6. The 14"x20" aperture was 
centered at a, 5 = 05 h 35 m 14!71, -05 o 23'4l"5) (equinox J2000). This position is 
the center of position 1SW, 18.5" S and 26.2" W of ^Ori C, observed with HST 
(see Rubin et al. 1997). The orientation of the long axis of the aperture was at 
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Figure 6. ISO SWS02 spectrum of the 2.8935 fim line in the Orion 
Nebula taken at position 1SW (see text for details). Observations were 
made with a 14"x20" aperture. The data (asterisks) have been fit with 
a Gaussian and a linear baseline. The residuals, data point minus curve 
fit, are indicated as the dashed line. 



PA 190.28°, determined by the time of the actual observation. With a correction 
for Yhelio = 17 km s _1 from emission lines arising in the main ionization zone 
(O'Dell et al. 1993), the line has a rest wavelength 2.89350±0.00003 fim (3456.02 
cm -1 ); it is a factor of 3.6 weaker than the nearby, H I 11-5 line (2.8728 /jm). 
The average surface brightness in this aperture was 3.9E-18 W m _2 arcsec -2 . 

We do not detect the uid line in our SWS02 spectra of any of the five bright 
PNs, which we observed for a comparable time. These are the PNs discussed 
earlier plus NGC 7009 (the Saturn Nebula). We have also examined many of the 
ISO archival spectra taken by others of the Orion Nebula as well as several other 
PNs. We find no indication of the uid line in the PNs examined, which included 
NGC 7027, NGC 6543, the low ionization PNs IC 418 and BD +30 3639, as well 
as the bright H n region M 17. Unlike our 1SW position (well in the heart of 
the ionized region), most of the other ISO Orion spectra were taken at positions 
where molecular/ photodissociation (PDR) /neutral species would be expected 
to dominate. There is some indication of the uid line in the following two archival 
spectra: 

• an SWS01 low-resolution scan of "Orion Bar d2" (PI Drapatz) a, 5 = 
05 h 35 m 21!40, -05°25'40'/l (J2000), taken 1997 October 11. The exposure was 
1912 sec. with a 20"x33" aperture. The long axis of the aperture was at PA = 
189.87°. We measure a line flux of 3.25E-19±0.75E-19 W cm" 2 . 

• an SWS06 scan of "Orion-BN" (PI Cernicharo) a, 5 = 05 h 35 m 14 f ;20, 
-05°22'23'.'6, taken 1997 October 12 (UT). The exposure was 7598 sec. with 
a 20"x27" aperture. The long axis of the aperture was at PA = 190.36°. We 
measure a line flux of 4.11E-20±1.22E-20 W cm" 2 . 
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Figure 7. Portion of UKIRT long-slit spectrum of the 2.8935 /im 
line in the Orion Nebula that is centered at position 1SW (see text for 
details). This 1-D plot represents the spectrum through a 0.74 arcsec 2 
aperture. Error bars are shown. 

Spectra of Orion taken at the United Kingdom Infrared Telescope (UKIRT) 
confirm the presence of this strong uid line. A long-slit observation (80" x 1") cen- 
tered at 1SW and pointing through ^ 1 Ori C shows that the emission is spatially 
extended and appears to be coincident with the brightest part of the H n region. 
The highest surface brightness occurs at 1SW: 1.2E-17 W m~ 2 arcsec~ 2 (±20%) 
with FWHM = 24 km s _1 (resolved). Fi gure 7 shows a spectrum at this location 
for which the wavelength scale has not been corrected for the Doppler velocity 
of the source. When the correction is made, the rest wavelength is 2.8935 /im. A 
second observation was made with the slit passing through the proplyd P159-350 
(O'Dell & Wen 1994) and pointing toward ^ 1 Ori C. This spectrum also shows 
extended emission in the uid line with no significant enhancement of the line at 
the proplyd location. Because P159-350 has N e at least as high as 10 5 cm~ 3 
as measured by the [C in] 1909/1907 flux ratio from GO-7514 HST STIS pro- 
gram (PI Rubin) (Henney & O'Dell 1999), it is a likely indication that the uid 
line emissivity is not strongly enhanced at higher electron density, which may 
indicate the line is collisionally excited and not a recombination line. 

We checked against lists of wavelengths of atomic ions and molecules. For 
this, the atomic line list provided at www.pa.uky.edu/~peter/atomic by Peter 
van Hoof was a particularly useful resource. There were two ionic lines that 
were very close to the uid wavelength and that were "reasonable" candidates. 
The first (best) candidate is: [Cr III] 2.89349±0.00020 ^m (3d 4 3 G 4 -» 3d 4 3 H 5 ) 
arising from 20851.87 cm -1 (~2.6 eV) above ground (Ekberg 1997). Chromium 
has a solar number abundance log [Cr/H] + 12 = 5.67 (Grevesse & Sauval 
1998). The LP. range for Cr++ is 16.49 - 30.96 eV, which would make Cr++ 
one of the dominant Cr species in Orion although the fraction of Cr +3 may be 
higher. However, we do not believe that this [Cr III] line is a likely identification 
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because there are other [Cr in] lines for which we have upper limits. These are 
[Cr m] 2.77830±0.00020 /an (3d 4 3 G 5 -» 3d 4 3 H 5 ) arising 20995.16 cm" 1 above 
ground; 2.88560±0.00022 (3d 4 3 G 5 -> 3d 4 3 # 6 ); and 2.91564±0.00020 fim 
(3d 4 — > 3d 4 3 #4) arising 20702.45 cm" 1 above ground. There was no hint 
of these lines in various SWS02 bandpasses that we observed covering these 3 
other transitions. From the similarity of the levels within the multiplet and A J 
of or -1, it is difficult to understand how the 2.89349 line would be so strong 
and all the others absent. 

The second close coincidence is an [Fe v] line at 2.8934±0.00075 /um (3d 3 £?5 
— > 3d 4 3 -F4) arising from 30430.10 cm" 1 above ground. Fe v is isoelectronic with 
Cr III. The LP. range for Fe +4 is 54.8 - 75.0 eV, which is higher than the He ++ 
limit of 54.4 eV. There is a stringent upper limit on the ionizing flux greater 
than 54.4 eV in Orion by virtue of deep spectroscopy that set a "record" limit 
on the flux ratio He n 4686/H/3 < 7xl0~ 5 (Baldwin et al. 2000). Thus, there 
should be negligible Fe v in Orion. 

Our assessment of what ionization range might produce the uid line is that 
it is likely from a species having a production LP. between 13.6 and 54.4 eV 
because it is not seen in the PNs we examined. However the uid line is also 
not detected in the H n region M 17 or the low-ionization PNs IC 418 and 
BD +30 3639. In Orion, the uid line is much stronger in the ionized region and 
only marginally present (if present at all) in the PDR/molecular /neutral zones. 
The fact that it is found in only the Orion Nebula reminds one of the uniqueness 
of proplyds to Orion (at least before they were found in some other nebulae). 

There is a reward offered for the correct identification of this uid line. The 
details were presented at the meeting and will be provided to anyone interested 
by writing to the first author. 



5. Upcoming and Planned IR Missions 

In this section, we discuss planned IR spectroscopic capability with a major em- 
phasis on NASA and ESA projects. There are some comparisons with the past 
projects/missions such as KAO, IRAS, and ISO. With one exception, we do not 
cover ground-based observatories but do remind the reader of excellent spectro- 
scopic capabilities to be afforded by these facilities (e.g., Gemini North/South, 
Keck, and the VLTs). 

5.1. Space Infrared Telescope Facility (SIRTF) 

This is a NASA space observatory with an 0.85-m telescope. It will have three 
cryogenically-cooled instruments. The launch date is expected in July 2002 
with a goal of a lifetime in excess of 5 years. SIRTF will be in an Earth- 
trailing, heliocentric orbit. The main spectroscopic instrument is the Infrared 
Spectrograph (IRS) (PI Jim Houck). There will be 4 separate modules that cover 
5.3 - 40 /um at resolution (R) from 60 - 120 (low-R mode) and a higher resolution 
(R = 600) capability from 10 - 37 /mi. It will provide a major advance in 
photometric sensitivity as can be seen in Figure 8 which shows comparisons with 
KAO, IRAS, ISO, and SOFIA. The angular resolution will be an improvement 
over ISO but somewhat less than the KAO as is seen in Figure 9. The URL 



http: / /sirtf.caltech.edu| provides further details. 
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5.2. Stratospheric Observatory for Infrared Astronomy (SOFIA) 

This is a joint project of NASA and the German aerospace center DLR. SOFIA 
will be the successor to the KAO. The facility is a modified Boeing 747SP air- 
plane with a 2.5-meter telescope. First light is expected in late 2004, with general 
observing commencing in 2005. The planned lifetime is 20 years and the home 
base is Moffett Field at the NASA/Ames Research Center. 

The photometric sensitivity is shown in Figure 8 in comparison with other 
missions. The angular resolution will be a substantial improvement, especially 
if the goal (solid curve) is attained (see Figure 9). The complement of First 
Light Instruments now being built are all listed below, but we describe here 
only those with medium to high spectral resolution capability. Figure 10 shows 
the instrument coverage, along with those of SIRTF, in R-A space. 
CAMERAS : There are 3 facility instruments - HAWC (PI Harper), FORCAST 
(PI Herter), and FLITECAM (PI McLean). FLITECAM will provide moderate 
resolution spectroscopy (R ~1000 - 2000) from 1 - 5.5 ^m. 
IMAGING PHOTOMETER : There is one Principal Investigator (PI) instru- 
ment, HIPO (PI Dunham). 

SPECTROGRAPHS : There are one facility instrument, (1) below, and five PI 
instruments (2)— (6) below - 

(1) the Airborne Infrared Echelle Spectrometer (AIRES) (PI Erickson) with 
ultimate design goals from 17 - 210 /xm with R ~ 10,000 (higher at 17 /im). 




Figure 8. The photometric, one-cr, point source sensitivities after one 
hour of integra tion for IRAS, KAO, ISO SIRTF, SOFIA, and Herschel. 
Graphics from http:/ /sofia. arc. nasa.gov. 
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Figure 9. The spatial resolutions (~50% encircled energy diameter) 
for IRAS, KAO, ISO, S IRTF, SOFIA, and Herschel. Graphics from 
http:/ /sofia. arc. nasa.gov. 



(2) the CAltech Submillimeter and far-Infrared Mixing Receiver (CASIMIR) 
(PI Zmuidzinas) that is a heterodyne receiver which will operate from 250 - 
600 /um with high spectral resolution. 

(3) the Echelon Cross Echelle Spectrometer (EXES) (PI Lacy) - a grating spec- 
trograph that will operate from 5-28 /im with high (R ~ 10 5 ), medium (R ~ 
10 4 ), and low (R ~ 3000) spectral resolution. 

(4) the Field Imaging Far-Infrared Line Spectrometer (FIFI LS) (PI Poglitsch). 
This will cover from 45 - 210 ^m with R ~ 1700. 

(5) the German Receiver for Astronomy at Terahertz Frequencies (GREAT) (PI 
Guesten). This is also a heterodyne receiver and will operate from 75 - 250 fim 
with high R. 

(6) the Submillimeter And Far InfraRed Experiment (S AFIRE) (PI Moseley). 
This is an imaging Fabry-Perot bolometer array spectrometer which will operate 
from 145 - 655 /xm with an ultimate goal of reaching R = 10,000. 

Further details are available via the SOFIA web site at 
http:/ /sofia. arc. nasa.gov/ 

5.3. Herschel Space Observatory (HERSCHEL) 

Formerly called FIRST, this is an ESA "Cornerstone Mission" with substantial 
NASA participation. It consists of a 3.5-m diameter telescope radiatively cooled 
to ~80 K with coverage in the 60 - 670 /mi range. After launch, expected to 
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Figure 10. The spectral coverage and resolution that will be possi- 
ble on SOFIA (first-light instruments only) and SIRTF. The science 
instruments shown ar e those for SOFIA, exce pt for IRS, IRAC, and 
MIPS. Graphics from [http:/ /sona,arc.nasa,gov| . 



be in 2007, it will be positioned to orbit around the Earth-Sun L2 Lagrangian 
point (located 1.5xl0 6 km in the anti-Sun direction from Earth). It is planned 
to operate for at least three years. There will be three instruments in the helium 
cooled focal plane: (1) SPIRE, a long- wavelength (200 - 670 /im) camera and 
low-resolution spectrometer; (2) PACS, a short-wavelength (60 - 210 /im) cam- 
era and low-resolution spectrometer (R~1700); and (3) HIFI, a high-resolution 
(heterodyne) spectrometer (R similar to CASIMIR), that will have contiguous 
coverage from 240 - 625 /mi (480 - 1250 GHz) plus a channel from 157 - 213 /im 
(1410 -1910 GHz). NASA is contributing detectors and components for two of 
the instruments and plans to support US science participation through a US 
Herschel Science Center at the Infrared Processing and Analysis Center (IPAC). 
For more detailed information see [http:// sci.esa.int/home/herscheT7| 



5.4. Planned IR Missions and the US Decadal Report 2001 2010 

Recently the latest US Astronomy Decadal Report was published entitled "As- 
tronomy and Astrophysics in the New Millennium" (National Research Coun- 
cil, NRC 2001). This often is referred to as the McKee/Taylor Report after 
the co-chairs of the Astronomy and Astrophysics Survey Committee. The Re- 
port can be accessed on-line at www.nap.edu/books/0309070317/html. Among 
other things, the report prioritized equipment initiatives (combining ground 
and space). Before continuing, we note that the committee reaffirmed the 
recommendations of the prior Survey Committee (NRC, 1991) and endorsed 
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the completion of SIRTF and SOFIA as well as other projects (NRC, 2001, p. 
24). 

The highest priority of the Major new initiatives is the Next Generation 
Space Telescope (NGST) (NRC 2001, p. 37). Below we briefly discuss NGST, 
as well as two other recommendations with major IR spectroscopic promise - 
the #2 and #7 ranked Major initiatives, the Giant Segmented Mirror Telescope 
(GSMT) and the Single Aperture Far Infrared (SAFIR) Observatory, respec- 
tively. 

NGST will be a segmented, filled-aperture 6-m class telescope with a large 
sunscreen. It will be in an orbit > 10 6 km from Earth to achieve very cold 
operating temperatures. The goal is for NGST to cover the range ~ 0.6 - 
27 fim. It is expected to be a factor of 10 4 - 100 more sensitive than SIRTF 
(from 4-20 fim). Spectral resolution is expected to be as high as 5000. 
GSMT will be a ground-based, 30-m class telescope. Through adaptive optics, 
it will have diffraction-limited imaging in atmospheric windows from 1-25 fim. 
The full wavelength goal is 0.3 - 25 p. Spectral resolution is expected to be 
as high as R = 10 5 . Unlike the space observatories planned, it will have the 
capability of adding new instruments. 

SAFIR will be an 8-m class, passively cooled space telescope that is planning 
to cover the spectrum from 30 - 300 /im. It is anticipated to provide 2 to 10 
times the spatial resolution of SIRTF, SOFIA, or Herschel. Even more dramatic 
will be SAFIR's increase in sensitivity compared with these upcoming missions, 
reaching several orders of magnitude depending on wavelength (see Fig. 3.2 on 
p. 103 of NRC 2001). Spectral resolution is expected to be as high as 1000. 

6. Summary and Conclusions 

With regard to our efforts to derive electron densities, there are several instances 
of the observed line flux ratio being clearly out of range of the theoretical pre- 
dictions using current atomic data. This is most dramatically illustrated with 
the [Ne v] (14.3/24.3) flux ratio for NGC 6818. These are both strong lines 
observed with signal-to-noise ratios better than 100. The observed flux ratio is 
0.707, which is significantly less than the predicted ratio using collision strengths 
from Lennon & Burke (1994). Because the line ratio near the low-iV e limit de- 
pends on collision strengths and not the transition probabilities (A- values) , this 
result may point toward a need to reexamine the collision strengths. 

A very rapid benefit of the type of interaction the Conference was to pro- 
mote occurred in this case. Two participants, Nigel Badnell and Keith Berring- 
ton, had made new computations for the Ne V effective collision strengths, 
including the relevant transitions. The former work had just been published 
(Griffin & Badnell 2000). Their paper provides the most detailed calculations 
to date. Their recommended values for each transition are the last entries in 
their Table 5. In the low-iV e limit, F(14.3)/F(24.3) should scale as (Ci 3 /(Ci 2 + 
C 13 ))x (24.3175/14.3217), where C 12 = T( 3 P - 3 A) exp(- X i2 / 'kT e ) and C 13 = 
T( 3 Po — 3 P2) exp(—xi3/kT e ). The xs are the energy level differences and k is the 
Boltzmann constant. For T e = 10,000 K, this gives a flux ratio 0.8090, while our 
detailed solution for the 5 lowest levels yields 0.8097. Although not material for 
this calculation, we did include T( 1 D2 — 1 So) (transition between lowest levels 
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4 and 5) that had been omitted in their publication. The missing entry is 0.545 
at 10,000 K (Badnell, private communication). We use the most accurate wave- 
lengths for the two [Ne v] IR lines measured from high spectral resolution ISO 
observations of a PN (Feuchtgruber et al. 1997) to provide energy levels 2 and 
3 in our code. Their paper is a fine example of how astrophysical measurements 
can actually do better than any laboratory measurement or quantal calculations 
to determine fundamental atomic data. 

The result of using the Griffin & Badnell (2000) Ts is gratifying, bringing 
closer agreement with the observations (see Table 6) although the NGC 6818 
observed ratio is still smaller than the low-iV e limit for T e 10,000 K or any 
reasonable T e (using adjacent 6,300 and 25,100 K entries in their table 5). 

We use CS's for S III, Ar III, and Ar v from Galavis et al. (1995); for Ne III 
and Mg v from Butler & Zeippen (1994). As a result of this Conference, we 
learned of the new calculations of Ne III CS's (McLaughlin & Bell 2000). When 
we use these instead of Butler k, Zeippen values, we find that the F(36.0) /F(15.5) 
ratio in the \ow-N e limit, for T e = 10,000 K, is increased slightly to 0.090. This 
makes the discrepancy with the observed ratios only slightly less (see Table 4). 

Recently Dinerstein (2001) identified two long-standing cases of IR uid 
lines - at 2.199 and 2.287 /an. These were seen in the spectra of PNs (e.g., 
Geballe, Burton, & Isaacman 1991 and references therein). She identified the 
former as [Kr in] 3 Pi- 3 P2, and the latter as [Se iv] 2 P^/2^ 2 Pi/2- The identifica- 
tion of such heavy elements in the IR suggests that there is much future work to 
be done in the Conference "fields of endeavor" . To extract astrophysical knowl- 
edge, e.g., nebular elemental abundances of such species, it will be necessary 
to observe/identify spectral lines, add to the atomic database for many more 
heavy elements, and add their treatment to photoionization codes. In addition, 
there can be expanded N e diagnostic capability by enlarging the number of ionic 
species observed that have 3 P ground states. 

As this work shows, progress in the future will also depend very much on 
the ability to obtain co-spatial observations of the line pairs in the IR used as an 
N e diagnostic. Thus it is important that future astronomical facilities be able to 
do this and maintain a level of flux calibration sufficient for analyses of electron 
density, such as we had hoped to perform. 
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